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ABSTRACT 

The 18806 ROSAT All Sky Survey Bright Source Catalog (RASS/BSC) X-ray 
sources are quantitatively cross-associated with near-infrared (NIR) sources from 
the Two Micron All Sky Survey Point Source Catalog (2MASS/PSC). An associ- 
ation catalog is presented, listing the most likely counterpart for each RASS/BSC 
source, the probability Pd that the NIR source and X-ray source are uniquely 
associated, and the probability P no -id that none of the 2MASS/PSC sources 
are associated with the X-ray source. The catalog includes 3853 high quality 
(P id > 0.98) X-ray-NIR matches, 2280 medium quality (0.98 > P id > 0.9) 
matches, and 4153 low quality (0.9 > Pid > 0.5) matches. Of the high qual- 
ity matches, 1418 are associations that are not listed in the SIMBAD database, 
and for which no high quality match with a USNO-A2 optical source was pre- 
sented for the RASS/BSC source in previous work. The present work offers a 
significant number of new associations with RASS/BSC objects that will require 
optical/NIR spectroscopy for classification. For example, of the 6133 Pd > 0.9 
2MASS/PSC counterparts presented in the association catalog, 2411 have no clas- 
sification listed in the SIMBAD database. These 2MASS/PSC sources will likely 
include scientifically useful examples of known source classes of X-ray emitters 
(white dwarfs, coronally active stars, active galactic nuclei), but may also con- 
tain previously unknown source classes. It is determined that all coronally active 
stars in the RASS/BSC should have a counterpart in the 2MASS/PSC, and that 
the unique association of these RASS/BSC sources with their NIR counterparts 
thus is confusion limited. 

Subject headings: catalogs - - infrared: general - - stars: coronae — X-rays: 
galaxies — X-rays: general — X-rays: stars 
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INTRODUCTION 



The ROSAT All Sky Survey Bright Source Catalog (RASS/BSC; IVoges et all Il999h 
contains flux and positional information for 18806 X-ray sources^, down to a limiting 0.1- 
2.4 keV count rate of 0.05 counts s _1 . The observations were made in 1990/91 with the 
ROSAT/PSPC satellite, and at a brightness limit of 0.1 counts s _1 (5843 sources) the 
RASS/BSC represents a sky coverage of 92%. The median positional uncertainty (la, in- 
cluding 6" systematic uncertainty) is 11". 

The sources in the RASS/BSC are the brightest X-ray sources in the sky. Brighter 
sources generally allow for better X-ray spectra and time variability observations, so the 
RASS/BSC sources are ideal for studying the X-ray properties of their respective classes 
with modern X-ray telescopes. However, as of 2008 Aug 1 only 1982 RASS/BSC sources 
have been observed with either Chandra or XMM-Newton. Further, around two thirds of 
the sources in the RASS/BSC remain unclassified (cf. §H), so the number of known source 
classes in the RASS/BSC and the number of sources of each class is limited. Classifying 
RASS/BSC sources therefore not only expands our knowledge of the brightest X-ray sources 
in the sky, but also provides ideal targets for follow-up observations that can improve our 
understanding of specific source classes. As a larger fraction of the RASS/BSC is classified, 
the catalog will also become a useful tool in producing population constraints for various 
classes of objects, and the sources that remain unclassified can be targeted in searches for 
previously unknown source classes and objects with unusual properties. 

To classify RASS/BSC sources, one needs spectral and flux measurements in other 
wavebands. The error cones of the X-ray sources often contain several off-band sources, any 
one (or none) of which may be a counterpart. If the probability of association between the 
X-ray source and one or more off-band counterparts is calculated, the classification of that 
X-ray-off-band pair /group can be constrained. Furthermore, probable off-band counterparts 
can be studied with ground-based or other instrumentation, which can provide constraints 
on the source of X-rays. Similarly, the absence of off-band sources near a RASS/BSC source 
can also constrain its class. Thus, searching for off-band counterparts to RASS/BSC sources 
is a good way to advance the classification and study of these sources. 

A summary of previous work identifying likely counterparts for > 500 RASS/BSC 
sources is given in Table [U which also includes the present work for comparison. Many 
of these previous cross-associations searched for counterparts in catalogs of objects of a 



1 The reference (jVoges et al Jl999h refers to 18811 sources, but 5 duplicates were removed in the most recent 



version of the catalog. See f tp : //f tp . xray .mpe .mpg .de/rosat/catalogues/rass-bsc/1 . 3_vs_l .4 . dif f 



for details (accessed 2009 May 5). 
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specific type, like stars (IHuensch et al.l Il998l; iHiinsch et al.l Il999l ; iMakarov Sz Urban! I200CH; 

Torres et al. 2006 ). galaxies ( Zimmermann et al.l 2001 ). or quasars ( Bade et al. 19981 ; Zickgraf et al 
20031 ). Other previous cross-associations searched for counterparts i n catalogs of sources 



in one or more spe c ific wave-bands, like optical (jVoges et all Il999l: 



Schwope et al.ll2000l ; iMcGlynn et al.l 120041 ) and/or radio (IBauer et al. 



Rutledee et a l. 2000; 



2000l : iMcGlvnn et al. 



2004 ). Only work that specifically addresses finding counterparts to RASS/BSC sources has 
been included in Table [TJ while work t hat includes other X-ray ca talogs without explicitly 
considering the RASS/BSC subset (e.g. iFlesch fc Hardcastlei 12004 ) has not been included. 



The RASS/BSC sou rces were quantitati vely cross-associated with USNO-A2 optical 
sources in previous work (IRutledge et al.ll2000l ). The present work adapts the methods used 
in that X-ray-optical cross-association to associate R ASS/BSC sources wit h Two Micron 
All Sky Survey Point Source Catalog (2MASS/PSC; ISkrutskie et al.l 120061 ) near-infrared 
(NIR) sources. The present work is the first to produce a cross-association between the full 
RASS/BSC and 2MASS/PSC catalogs. 

2MASS surveyed 99.998% of the celestial sphere in the near-infrared, resulting in a Point 
Source Catalog containing 470992970 point sources, and an Extended Source Catalog (XSC) 
containing 1647599 extended sources. The majority of 2MASS/XSC sources are galaxies 
(-97%), while most 2MASS/PSC sources are stars in the Milky Way. The 2MASS/PSC 
also contains point-source processed versions of virt ually all the 2MA SS/XSC sources, as 
well as a significant number of unresolved galaxies (ICutri et al.ll2003l ). The present work 
considers only the 2MASS/PSC, for which the greatest sensitivity was ac hieved in the J 
band (1.12 — 1.36/mi), with a 10cr detection level better than 15.8 mag (ISkrutskie et al. 



2006|). 



This paper is organized as follows. In §2J the method used to perform the statistical 
cross- association between RASS/BSC and 2MASS/PSC sources is described. The results of 
this cross-association are discussed in £j3l which includes an association catalog containing 
the most likely NIR counterpart for each X-ray source. In §H the source types listed in the 
SIMBAD database for the 2MASS/PSC sources that appear in the association catalog are 
discussed. The present work is compared to previous cross- associations in §5J and discussion 
and conclusions can be found in 3B1 



2. ANALYSIS 



The present approach follows that of previous work, in which RASS/BSC X-r ay sources 
were statistically associated with USNO-A2 optical sources (IRutledge et al.ll2000i . R00 here- 
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after). The main steps to this approach are recounted here, as relevant to the present work, 
while the earlier work is referred to for a more detailed discussion. 

Source fields are defined as regions within 75" of a RASS/BSC X-ray source. For each 
source field, a number of non-overlapping, adjacent background fields are extracted (in R00, 
10 fields were extracted along a line of constant declination through the RASS/BSC position). 

One defines a function, LR, which depends on the joint properties of an X-ray-off-band 
source pair, such as the X-ray-off-band source separation r and its uncertainty, and the 
off-band magnitude of the candidate counterpart. LR should be large when the the X-ray 
source and off-band counterpart are associated, and small when they are not. 

LR is computed for each X-ray-off-band source pair within a source field, and for off- 
band sources in background fields as if there is an X-ray source at the background field 
center. 

The reliability of association Ri of an X-ray source with a given off-band source is 

R{LRi) = n trne(LRi) 

ThTue(LRi) + nf a i se (L.Rj) 

where LRi is the value of LR for that X-ray-off-band source pair, n tTne (LRi) is the number 
of true associations (per source field) with that LR value, and n{ a \ se (LRi) is the number 
of associations (per source field) with that LR value, where the off-band source is an un- 
related background source. As it is not known a priori which off-band sources in source 
fields are unrelated background objects and which are true counterparts to a RASS/BSC 
source, R(LR) must be approximated based on the number distributions in LR of off-band 
sources per source and background field, n sonrce (LR) and ^background (LR) respectively. The 
approximation to R(LR) used is (cf. R00) 

7->/ r 7-,\ ""source 

(LR) — "^background (-^-R) / Q x 
R(LR) W — — . (2) 

"■source 

In the theoretical development below, the formal definition of R (eq. [Tj) is used, while the 
approximation to R (eq. [2]) is used whenever discussing computed values or results. 

Assuming that each X-ray source is either associated uniquely with one of the M off- 
band sources in the corresponding source field or associated with none of them, the proba- 
bility of unique (exclusive) association between the X-ray source and the candidate off-band 
counterpart, Pid,i, is calculated: 



(3) 
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Further, the probability that none of the off-band sources are associated with a given X-ray 
source, 

-Pno-id = R d ' (4) 

1 + 2^j=l 1-Rj 

is computed. This permits the fraction of the RASS/BSC that has an identifiable counterpart 
in the off-band catalog (the 'quality') to be determined: 

Q=(l- Pno-id) , (5) 

where the average is over all RASS/BSC sources. The definitions of and P no -id given 
above are the same as in the previous work (R00). The apparent difference is due to the 
normalization factor being written out explicitly in the present work, with common factors 
being cancelled from the numerator and denominator. 

The remainder of §0 describes the adaptation of the method described above to the 
RASS/BSC-2MASS/PSC cross-association performed in the present work. The main focus 
of this analysis is to identify individual 2MASS/PSC sources associated with a RASS/BSC 
source. 



2.1. Source and Background Fields 

NIR sources are extracted from the 2MASS/PSC catalog for source and background 
fields, within a radius of 75" from each field center. The background fields are centered on a 
3x3 square lattice of spacing 150" around each source field, and assigned the same positional 
uncertainty as the X-ray source in the source field. The square lattice is chosen to ensure 
that the background fields are taken from regions similar to the source fields. Background 
fields containing a RASS/BSC source are excluded from the analysis, since it is assumed in 
calculating R(LR) that background fields contain no NIR sources that are associated with 
a RASS/BSC source. There are a total of 150132 background fields. 

The distance from a 2MASS/PSC source to the field center is denoted r. The differential 
number distributions of separations r in the source and background fields are shown in Fig. [TJ 
At each r, the shapes of the cumulative distributions from r to 75" are compared using a 



two sample Kolmogorov-Smirnov (K-S; iPress et al.l Il992l ) test. The K-S test gives the 



probability that the cumulative source and background field distributions arise from the 
same underlying distribution. For r < 30" there is an excess of NIR sources in the source 
fields compared to the background fields, and the K-S probability is below 10~ 3 for r < 25". 
This excess in the source fields is due to 2MASS/PSC sources that are associated with a 
RASS/BSC source. 
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There is an improbable difference (p ~ 4-10 -4 ) between the number distributions of NIR 
sources, observed when comparing source and background fields in the range 71" < r < 75". 
A detailed examination of these distributions finds a small, but apparently significant, excess 
of NIR sources in the source fields compared to background fields. Only 314 RASS/BSC 
sources have positional uncertainties greater than 30", so it is unlikely that the excess of 
NIR sources at r ~ 71" is due to 2MASS/PSC sources associated with a RASS/BSC source. 
However, the clustering of 2MASS/PSC sources (in open clusters or galaxy clusters) could 
cause deviations in the source field distribution compared to that of background fields, for 
r > 30", associated with the presence of the X-ray source. 

Differential number distributions in J-band magnitude wij for NIR sources in the source 
and background fields are shown in Fig. [2J There is an excess of NIR sources in the source 
fields, except for mj > 17.6, where there is a slight deficit of NIR sources in the source fields 
compared to the background fields. This deficit is likely due to the detectability of faint NIR 
sources being limited by the presence of a larger number of bright NIR sources in the source 
fields. The excess of NIR sources in source fields is statistically most significant for bright 
(mj < 10.6) NIR sources. 

The excess NIR sources present in source fields compared to background fields are NIR 
sources that are associated in some way with the X-ray source. A large number of known 
X-ray source classes are associated with single, point-source NIR counterparts, like coronally 
active stars or (spectroscopic or eclipsing) binaries, active galactic nuclei (AGN), or low-mass 
X-ray binaries. The present analysis is optimized to find these sources classes, as they have 
single NIR counterparts. However, there are known X-ray source classes which are associated 
with multiple NIR sources, either due to the nature of the source class (e.g. X-ray clusters of 
galaxies) or because the source class tends to be found in groups of similar NIR sources (e.g. 
T Tauri stars in stellar open clusters, or AGN in galaxy clusters). Further, X-ray emitting 
NIR resolved binaries also have more than one NIR source associated with the X-ray source. 
These classes of X-ray sources that are associated with multiple NIR sources will still appear 
in the present analysis, but the top Pd match will be affected by source confusion; multiple 
possible NIR counterparts will decrease the Pid value for the most likely counterpart. The 
remainder of this section attempts to constrain the number of RASS/BSC sources that are 
associated with a given number of NIR sources. 

The normalized number distribution of RASS/BSC sources with Anir NIR sources 
in the corresponding source field is shown in Fig. [3j for Anir < 50. The distribution of 
RASS/BSC sources has a sharp peak around a mode of A NI r = 3, and has a tail extending 
to Anir = 373. There are 1199 RASS/BSC sources with 50 < Anir < 373, with less than 
30 RASS/BSC sources at any given Anir. 
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A source field containing a RASS/BSC source that is associated with m NIR sources 
should on average contain m more NIR sources than an adjacent background field. Thus, the 
normalized distribution in iV NIR of RASS/BSC sources associated with m NIR sources should 
be consistent with the normalized number distribution of background fields in N^m — m. 
Further, the overall distribution of RASS/BSC sources in A^ir should be consistent with a 
weighted sum over m of background field distributions in A^nir — m, where the weights are the 
fraction of RASS/BSC sources that are associated with m NIR source s. To estimate th ese 
weights, a Levenberg-Marquardt nonlinear x 2 minimization fitting (cf. Press et alii 19921 ) to 



the Anir < 50 distribution of RASS/BSC sources is performed using the Origin (OriginLab, 
Northampton, MA) software. Values of iV NIR > 50 are excluded since the uncertainties are 
assumed to be Gaussian, which is not necessarily a good approximation for the low number 
counts in this range of Aj^r. Since there are almost eight times as many background fields as 
RASS/BSC sources, the fractional uncertainties on the background field distribution should 
in general be less than half those on the RASS/BSC distribution, though differences in the 
distributions may increase the importance of the background field distribution uncertainties. 
The background field distribution uncertainties are neglected when fitting, and it can be 
verified after fitting that these uncertainties are indeed negligible. 

The aim of the fitting is to determine whether there is an acceptable description of the 
observed distribution of RASS/BSC sources in A^ir, under the assumption that a fraction 
f m of all RASS/BSC sources have m NIR counterparts, where m = 1, 2, 3, ... . The method 
adopted is to first fix f m = for all m except m — 0, and then free f m components one by 
one (first m — 1, then m — 2, etc.) until the reduced x 2 drops below 1. For the f m that 
are allowed to vary in each step, a constraint f m > 10~ 10 is adopted to avoid computational 
problems that arise in the fitting routine if the constraint f m > is used. With the addition 
of the m = 7 component, the reduced x 2 value is 0.94, so no further components are added. 
The m = 4, 5, &6 components have best-fit values f m = 10~ 10 , i.e. equal to the lower 
constraint, so these are then fixed to zero. Re- running the fit, the m = 3 component is 
consistent with zero weight (/3 = 0.02 ± 0.04), so this weight is also fixed to zero. The final 
best fit (reduced x 2 = 0.86) is then /„ = 0.12 ± 0.01, ft = 0.66 ± 0.03, f 2 = 0.18 ± 0.03, 
and ft = 0.040 ± 0.006. Thus, 78 ± 3% of all RASS/BSC sources have zero or one NIR 
counterparts, i.e. satisfy the assumptions used when calculating P id and P no -id- Extending 
the present method to include matches with two or more NIR sources in the calculation of 
Pid will be the focus of future work. 

It should be noted that the quoted uncertainties on f m only take into account the 
diagonal elements of the covariance matrix, and thus do not consider dependencies between 
different m components. Since the different components only differ by a shift in A^nir of 
1, f m parameters with similar m are likely to be highly dependent, so the true fraction of 
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the RASS/BSC that is associated with m NIR counterparts may not be consistent with the 
quoted uncertainty on f m . Further, it is likely that the m = 7 component accounts for a 
population of RASS/BSC sources that are associated with 5 < m < 9 NIR sources, rather 
than just an m = 7 population. 

Given the best-fit values of f m , it can now be verified that neglecting the uncertainties 
on the background field distribution in iV NIR is an appropriate assumption. The largest ratio 
of the uncertainty on the RASS/BSC distribution to that on any of the f m components is 
0.27 (for m = 1 and iV NIR = 34). Adding this uncertainty to the RASS/BSC distribution 
uncertainty in quadrature would give an uncertainty that is 4% larger than the RASS/BSC 
distribution uncertainty alone, so neglecting the background field distribution uncertainty is 
unlikely to have a significant effect on the fitting. 



2.2. LR 

Large uncertainties in the RASS/BSC positions make it difficult to identify likely NIR 
counterparts based only on the probability of positional coincidence. Therefore, the J-band 
brightness of the NIR source and the local density of NIR sources were included in the 
function LR: 

_ exp(gf)exp(-p7rrf) 
LKl ~ 2vra 2 iV(< mjli ) • { ) 

N(< mj t i) is the number of brighter NIR sources in background fields, a 2 is the sum of 
squares of the uncertainties in the X-ray and NIR positions, and p is the the local density 
of NIR sources. To compute p, the number of 2MASS/PSC sources within 30' of the X-ray 
source is divided by angular area of a 30' radius cone. Some X-ray sources have a listed 
positional uncertainty of zero, for which an uncertainty of 12" is adopted. 

The quantity LR defined in equation (jBJ) should be understood as a figu re of merit, and 



not as a likelihood ratio in the formal statistical definition of that term (cf. IWilkslll938l ): in 
the formal definition the likelihood ratio is a ratio of probability density functions, whereas 
the quantity LR as used here and previously (R00) includes integrated probability density 
functions, and thus does not satisfy this definition. Moreover, the quantity LR is not here 
used as a likelihood ratio, in that a likelihood ratio is used to determine the relative likelihood 
of two competing hypotheses; here, LR is used strictly as a figure of merit and not, itself, to 
compare the relative likelihood of two competing hypotheses 



1 



2 The quantity LR evolved (R00) as a modified version of a quantity referred to as a "likelihood ratio" , 
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The LR used in the present work differs from that employed in the RASS/BSC- 
USNO-A2 cross-association by a factor of (27rer) _1 , due to the use of a two dimensional — 
rather than one dimensional — normalization of the probability of coincidence, and by the 
added factor of exp(— pirrf), which is the probability that there is no closer NIR source to 
the field center, p ranges from 0.3 to 50 NIR sources per square arcmin, is 2.8 sources per 
square arcmin on average, and is typically large close to the galactic plane and small away 
from it. Including a p dependence in LR accounts for the fact that the higher the local 
density of NIR sources, the larger the chance that an unrelated background NIR source falls 
within the error cone of an X-ray source. 

Fig. H] shows the number distributions in LR of NIR sources in source and background 
fields, and the resulting R(LR). To ensure sufficient counting statistics, the number distri- 
butions are calculated by first dividing the \n(LR) range into bins of width 0.01, and then 
for each bin averaging over a wide enough range of bins around it — at least 21 — to keep 
the propagated Gaussian counting uncertainty ctr < 0.01, and to keep the total number of 
NIR sources above 250 — to justify using the Gaussian uncertainty. The normalization of 
the number distributions is such that J(N/&eld)d(ha.(LR)) approximately equals the total 
number of NIR sources per field for that field type (the equality is not exact because of 
the averaging over bins described above). R(LR) is set to zero for LR values less than the 
largest LR for which R(LR) = 0, so that R > for all X-ray-NIR pairs. Negative R values 
can arise from random fluctuations, or from breakdown of the assumption that source fields 
contain both background NIR sources and NIR sources that are associated with the X-ray 
source. An example of such breakdown is the apparent deficit of very faint NIR sources in 
source fields compared to background fields. 

X-ray-NIR pairs with In(LRi) values in the same bin are all assigned the same Ri 
value, so Ri values in the same ln(LRj) bin are perfectly correlated. Further, the averaging 
over bins used to calculate R(LR) introduces a correlation between bins. All bins with 
ln(LR) > —16.05 are correlated to some extent, since there are only 247 NIR sources in 
background fields with LR in this range, and any bin in that range of LR thus must include 
at least one bin with ln(LR) < —16.05 in the averaging described above. For the remaining 



used in previous work ( Wolstencroft et al. 19861 ) to identify likely off -band counterparts of IRAS sources. 
The terminology "likelihood ratio" was applied in the previous work | Wolstencroft et al.|[l986h . where the 
relevant quantity was in fact used to choose between two competing hypotheses; however, the quantity to 
which it was applied involved terms which were not probability density functions, so the quantity used there 
does not satisfy the formal definition of a likelihood ratio. Thus, the reader is strongly cautioned that in the 
present work, the quantity LR is used strictly as a figure of merit; it is not itself used as a likelihood ratio to 
distinguish between two competing hypotheses, and should not be considered as anything other than strictly 
a figure of merit. 
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range of LR where R(LR) is non-zero, —34.07 < ln(LR) < —16.05, the distance in ln(LR) 
over which a correlation exists varies between 0.23 and 2.33. There are 10191 source fields 
with at least one pair of matches for which the R values are not independent. Thus, care must 
be taken when propagating or to calculate opy, and if the uncertainties are propagated as if 
the R values are independent, it must be verified that this assumption does not significantly 
affect the crpid values of interest. 

The chosen form of LR neglects any dependence on p of the number distribution of 
NIR sources as a function of toj. Such a dependence could arise from source population 
differences between the denser galactic plane and the less dense regions at higher galactic 
latitudes, or from the detectability of faint sources being limited by the presence of bright 
sources. Some dependence on p is seen in a normalized version of iV(< mj ti ), and this is 
most significant for the very brightest sources (raj < 4). However, there are very few mj < 4 
sources in source or background fields (cf. Fig. [2]), and these are likely to have R values close 
to 1 when they appear in source fields (cf. §3§, so it is not expected that addressing these 
population differences would have a significant impact on the resulting association catalog. 



2.3. Control Fields 

In addition to the source and background fields, one control field (75" radius) is extracted 
per square degree of the celestial sphere. Control fields that overlap with a source field are 
excluded, bringing the total number of control fields to 41192. Each control field is assigned 
the properties of a randomly chosen RASS/BSC source, and LR i: R(LRi), and P^i are 
calculated for each NIR source in a given control field using the R(LR) function computed 
from the source and background fields. Note that R(LRi) is a function which depends only 
on the source and background fields, not on the control fields. 

A control field resembles an X-ray source that is not associated with a 2MASS/PSC 
source, so the P no -id values of control fields should be comparable to those of near-infrared 
faint RASS/BSC sources, and to RASS/BSC sources associated with any extended NIR 
sources that do not appear in the 2MASS/PSC. This information has prove n useful in pop 



ulation studies of blank- field X-ray sources, such as isolated neutron stars (IRutledge et al. 



20031 : iTurner et all 120091 ) . 
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3. RESULTS 

For 100 RASS/BSC sources, there are no 2MASS/PSC sources within 75". For the 
remaining 18706 X-ray sources, there are a total of 287736 distinct NIR sources within 
75". Except for three NIR sources that lie within 75" of both 1RXS J020220.5-010609 and 
1RXS J020221. 0-01071, each NIR source is paired with only one RASS/BSC source, so 
the distinction between an NIR source and an X-ray-NIR pair can be relaxed. Assuming 
Gaussian counting uncertainty on the number of NIR sources, there are 15.30 ± 0.03 NIR 
sources per source field, and 14.18 ± 0.01 per background field, amounting to an excess in 
source fields of 1.12 ± 0.03 NIR sources per field. 

The distribution in LR of the excess of NIR sources in source fields can be seen from 
Fig. HI For hi(LR) < —34.07, the source and background field distributions are consistent 
with no excess in source fields (0.03 ± 0.03 sources per field), and the probability that the 
same underlying distribution would give rise to two distributions as or more different than 
the sourc e and background field distributions — for \ia(LR) < —34.07 — is 0.55 (two sample 



K-S test; iPress et al.lll992l ). Thus, the number distributions of NIR sources in source and 
background fields are consistent with being identical for \n(LR) < —34.07, and setting 
R(LR) = for this range of LR is justified. For \n(LR) > -34.07, the excess of NIR 
sources in source fields compared to background fields is 1.10 ± 0.02 per field. Thus, there 
are 20600 ± 300 NIR sources in source fields in excess of the number of expected background 
sources for ln(LR) > —34.07, and these NIR sources are assumed to be associated with the 
RASS/BSC sources when calculating R(LR). 

Background fields are positioned near source fields to ensure that the population of 
NIR sources in background fields is similar to the population of background NIR sources in 
source fields, since the two populations are assumed to be identical when calculating R(LR). 
Further, the positional uncertainty of the RASS/BSC source in a given source field is adopted 
for the adjacent background fields, so that any anisotropy in the distribution of RASS/BSC 
positional uncertainties also is present for the background fields. The demonstrated similarity 
of the source and background field distributions of NIR sources with ln(L_R) < —34.07 
indicates that the background field NIR population does indeed approximate the background 
population in source fields, most likely also for ln(LR) > —34.07. 

There are 15.62 ± 0.02 NIR sources per control field, so there are more NIR sources per 
control field than per source field or per background field. Table Ogives the number of source 
and control fields in six different regions of the sky, and the average value of p — the local 
density of NIR sources — for the fields in that region. The regions are chosen to differentiate 
between the galactic poles, center, and anti-center, since these have different average NIR 
source densities in the 2MASS/PSC catalog. A larger fraction of control fields than source 
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fields lie in the galactic plane, where the NIR source density is larger than in the remaining 
part of the sky, and this explains why there are more NIR sources per control field than per 
source field or per background field. Only NIR sources with ln(LP) > —34.07 are of interest 
in the case of control fields, since only these sources have R > 0, and for this range of LR 
the cumulative distributions of NIR sources in control and background fields are consistent 
with arising from the same underlying distribution (two sample K-S test probability of 0.19; 



Press et al.lll992l ). The distributions thus only differ in normalization for this range of LR. 



To examine the impact of using Pd as a measure of the significance of an association 
rather than R, it is useful to consider the differential number distributions in R and Pd of 
NIR sources in the source and control fields. The distributions in R are shown in Fig. (5J 
while the distributions in pd are shown in Fig. [HI For R > 0.10 the number density in source 
fields is above that in control fields, while for R < 0.10 the number density in source fields is 
below that in control fields in all but two bins. In total, there is a 9% excess of NIR sources 
in control fields compared to source fields for R < 0.10, which is due to the larger average 
number of NIR sources in control fields. For R > 0.10, there is a 78% excess of NIR sources 
in source fields compared to control fields, and this excess is due to NIR sources that are 
associated with a RASS/BSC source only being present in source fields. The distributions 
in pd show that for Py > 0.28 the number density in source fields is above that in control 
fields, while for Pd < 0.28 the number density in source fields is below that in control fields 
in all but two bins. In total, there is a 6% excess of NIR sources in control fields compared 
to source fields for Pd < 0.28, and a 222% excess in source fields compared to control fields 
for pd > 0.28. If each field contained only one R > source, Pd would equal R for all 
NIR sources, and the distributions in pd would be identical to the distributions in R (cf. 
eq. [3]). However, 85% of source fields and 58% of control fields contain more than one R > 
source, so the distributions do in fact differ. The reason that a larger fraction of source fields 
than control fields contain multiple NIR sources with R > is because source fields contain 
NIR sources that are associated with a RASS/BSC source, while control fields do not. As 
a result, unrelated background sources in source fields are assigned lower pd values than 
similar sources in control fields, and it is possible that n contro i(Pd) > "-source (Pd) for some 
Pd = R where n contro i(P) < n source (P), as is the case for 0.10 < Pd < 0.28. 

The quality of the association catalog (eq. [5]) is Q = 0.70 ± 0.02. The uncertainty 
on Q is a conservative estimate of the statistical uncertainty, computed by noting that Q 
is simply the sum of of all Pd values in source fields divided by the number of RASS/BSC 
sources. Since each P id value is correlated with at least 249 other source field Pd values (due 
to the binning described in §2.2p . their erpid values cannot simply be summed in quadrature 
as if the Pd values are independent. Rather, the conservative assumption that all Pd values 
are perfectly correlated is adopted, and the uncertainties summed absolutely to obtain the 
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conservative estimate of the statistical uncertainty on Q. 

The interpretation of Q is that it is the fraction of the RASS/BSC with a counterpart 
in the 2MASS/PSC identifiable by the present method (R00). Therefore, 1 — Q is the 
fraction of the RASS/BSC that does not have an NIR counterpart, under the assumption 
that all RASS/BSC sources have zero or one NIR counterparts. However, it was shown in 
§2.11 that this assumption only holds for 78 ± 3% of RASS/BSC sources, so the assumption 
may introduce a large systematic uncertainty on 1 — Q. In fact, the fraction of RASS/BSC 
sources found in §2.1l to have no NIR counterpart is 12 ±1%, which is not consistent with the 
statistical uncertainty on 1 — Q. Thus Q, and by extention many P no -id values, have large 
systematic uncertainties introduced by the assumption of zero or one NIR counterparts to 
each RASS/BSC source. These systematic uncertainties are also likely to be present in the 
previous work (R00), as it was concluded there that ~ 27.2% of RASS/BSC sources do not 
satisfy the assumption that each X-ray source has zero or one USNO-A2 counterparts. 

Control fields do not contain NIR sources that are associated with a RASS/BSC source, 
but have P id and P no -id values computed based on the R(LR) function computed using the 
source and background fields. The control fields are thus like the source fields of X-ray bright 
NIR faint sources in the RASS/BSC. Fig. [7] shows the normalized cumulative distributions 
in 1 — P n o-id of RASS/BSC sources and control fields, n(< 1 — P no -id), for the six different 
regions of the sky described in TableEJ At every value of 1 — P no -id, ""-control fioids(< 1 — Pno-id) 
is larger for each region than for any other region with larger average NIR source density. 
Thus, X-ray bright NIR faint sources are more easily identified in regions of lower NIR source 
density. In the case of RASS/BSC sources, the distributions for regions in the galactic plane 
follow each other closely for 1 — P no -id < 0.75, and lie entirely below the distributions for 
regions outside the galactic plane. This may in part be due to greater confusion in the 
galactic plane as a result of higher NIR source density, but is most likely indicative of a 
difference in the population of X-ray sources in and outside the galactic plane. 



3.1. The Association Catalog 

Of the 18806 RASS/BSC sources, 18568 are associated with at least one 2MASS/PSC 
source with P id > 0. Further, each RASS/BSC source is assigned a P no -id- The association 
catalog (Table [3]) contains all RASS/BSC sources, their P no -id, and the 2MASS/PSC source 
in the corresponding source field with the highest non-zero Pd- If any of the RASS/BSC 
sources or 2MASS/PSC sources in the association catalog were listed in the SIMBAD 
database as of 2009 April 23, the source type given there is also included (cf. §U for de- 
tails). The association catalog is sorted by the 1RXS designation of the RASS/BSC source. 
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An online version of the catalog is also available!, and is searchable by p d and P no -id- 

The catalog contains 3853 2MASS sources with P id > 0.98, 2280 with 0.98 > P id > 0.9, 
and 4153 with 0.9 > p d > 0.5. These subsets of the catalog are referred to as high, 
medium, and low quality matches, respectively, and correspond to the catalogs which resulted 
from the RASS/BSC-USNO-A2 cross- association (R00). The background contamination 
(iV bk g = N — SjP idi , cf. R00) is 39, 101, and 1443 for the high, medium, and low quality 
matches respectively. 

The distributions of statistical uncertainties <7pi d on p d are shown in Fig. [HJ The high 
and medium quality matches all have opi d < 0.01, while the low quality matches have o~Pi d 
values clustered around 0.01, with a tail in the distribution extending to 0.03. For P id < 0.5 
almost all matches have crpi d values clustered around 0.01. The clustering around opi d ~ 0.01 
is due to the bin size used in calculating R(LR) being adjusted until <tr < 0.01, since the 
uncertainty in R is approximately the uncertainty in p d for sources where R ~ p d . The tail 
of larger opi d consists of sources for which there is a significant difference between R and p d 
due to the presence of other sources with non-zero R values in the same field. 

In propagating the uncertainty or to opi d , it is assumed that the R values in that field 
are independent. However, due to the method used to calculate R(LR), that may not be the 
case (cf. §2.ip . Thus, it must be verified that treating the R values as independent has not 
led to significantly smaller op; d values than those that would result from an error analysis 
considering correlations. An investigation of the matches in the association catalog reveals 
that none of the high or medium quality matches have R values that are correlated with 
another R value in the field, while 445 low quality matches and 2799 p d < 0.5 matches do 
have R values that to some extent are correlated with at least one other R value in the field. 
All of these low quality matches have R values that are correlated with exactly one other 
R value in the field, while 262 of the p d < 0.5 matches have R values that are correlated 
with between two and five other R values. A conservative rough estimate of the effect of 
correlations in the error propagation is that 0pi d is a weighted absolute sum rather than the 
square root of a weighted sum of squares, so considering correlations could increase op id by 
as much as a factor of square root the number of correlated R values. Thus, the uncertainties 
shown for the high and medium quality matches in Fig. [8] should be largely unaffected by the 
assumption that the R values are independent, while 11% of the low quality matches may 
have quoted uncertainties smaller by as much as a factor of 1/a/2 than those that would 
result from a covariant error analysis. Further, 31% of the p d < 0.5 matches may have 
uncertainties that are underestimated by as much as a factor of l/y/2, while 3% may be 



3 http : //dualcore . physics . mcgill . ca/RASS/ index . html 
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underestimated by as much as a factor of 1/ \^6. Based on the distributions shown in Fig. [HJ 
almost all low quality and P;d < 0.5 matches would thus have crpid < 0.03, also in a covariant 
error analysis. Whether the uncertainty on a low quality or Pd < 0.5 match is 0.01 or 0.02 
is unlikely to affect a decision to do follow-up observations of the NIR or X-ray source, nor 
is it likely to affect any conclusions drawn on the basis of the association, so it is concluded 
that the assumption of independent R values is adequate for the present work. 

In principle, correlations between two other R values in the field of a high or medium 
quality match could lead to a larger apid value in an error analysis considering correlations. 
However, those R values are likely to be small compared to the R value of the association 
catalog match (otherwise it would not be a high or medium quality match), so the contri- 
bution of the uncertainties on the other R values to crpid is likely to be small, whether or 
not the other R values are correlated. Thus, the conclusion remains that the uncertainties 
shown for the high and medium quality matches in Fig. [HJ should be largely unaffected by 
the assumption that the R values are independent. However, the uncertainties for the low 
quality and Pd < 0.5 matches may be affected by correlations between other R values in 
the field, so the fractions of these matches that are affected may be larger than those quoted 
above, as may the size of the effect in some cases. Fields for which there are many correlated 
R values are likely to be fields for which the assumption of one or no NIR counterparts to 
the RASS/BSC source is not valid, so the systematic uncertainty introduced by this assump- 
tion is likely to be larger than the propagated statistical uncertainty for these fields. Thus, 
whether the uncertainty is propagated as if the variables are independent or not would not 
significantly affect the quantitative or qualitative conclusions under the present analysis. 

The distribution of separations r between the X-ray and NIR sources in the association 
catalog is shown in Fig. [91 The distribution for high quality matches is the most sharply 
peaked, and a broadening is seen as Pd decreases. High, medium, and low quality matches 
are identified out to a distance comparable to the positional uncertainty for the RASS/BSC 
sources; the mode and average RASS/BSC positional uncertainty are 10" and 12.7" respec- 
tively. Several Pd < 0.5 matches are identified well beyond this range, some at distances 
greater than 45". 

Fig. [HI] shows the distribution of J-band magnitudes among the matches in the catalog. 
Almost all high and medium quality matches have mj < 12, while the low quality and 
remaining matches include a population of significantly fainter sources. As will be discussed 
in §4.11 many of the mj < 12 sources are coronally active stars, and the sharp drop in the 
number of matches at mj ~ 10 is likely due to coronally active stars fainter than this mj 
having an X-ray flux below the RASS/BSC detection limit. The second, J-band fainter, 
population of matches seen in Fig. [10] appears to match the brightness distribution of NIR 
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sources in the background fields, indicating that there is no strong X-ray to J-band flux 
correlation for these matches. 



4. SOURCE CLASSES IN SIMBAD 

The SIMBAD databas^l lists source types for many astronomical objects, including — as 
of 2009 April 23— all 18806 RASS/BSC sources, and 6115 of the 2MASS/PSC sources in the 
association catalog (Table [3]). SIMBAD is an incomplete, heterogeneous database, but offers 
a sense of the extent to which the association catalog includes previously unknown informa- 
tion. Objects listed in SIMBAD with the source type 'X-ray source' or 'Infra-Red source' are 
considered unclassified, and these include 12213 (65%) of the RASS/BSC sources and 229 of 
the aforementioned 2MASS/PSC sources. The remaining 6553 RASS/BSC sources and 5886 
2MASS/PSC sources are considered classified in SIMBAD. Of the 18568 matches in the as- 
sociation catalog, 2609 appear as associated in SIMBAD (the RASS/BSC and 2MASS/PSC 
identifiers are both listed for the same object). This section examines the source types and 
RASS/BSC-2MASS/PSC associations listed in SIMBAD for the 2MASS/PSC sources in the 
association catalog. 

All RASS/BSC-2MASS/PSC matches fall into one of four mutually exclusive categories: 
a RASS/BSC source and a 2MASS/PSC source may or may not be associated in SIMBAD, 
and the 2MASS/PSC source may or may not be classified in SIMBAD. The number of 
matches in the association catalog that fall into each of these categories is given in Table HJ 
for several different Pid cutoffs. For the classified 2MASS/PSC sources that are not associated 
with the corresponding RASS/BSC source in SIMBAD, the 2MASS/PSC source type can be 
adopted for the RASS/BSC source with confidence P;d- Moreover, this table points out where 
opportunities to classify RASS/BSC sources exist, through observations of the 2MASS/PSC 
source. For example, there are 1203 high quality matches for which the 2MASS/PSC source 
is not classified in SIMBAD, and follow-up optical/IR spectroscopy of these 2MASS/PSC 
sources could classify the X-ray source as an AGN, star, white dwarf, or other class of X-ray 
source. 

The number of 2MASS/PSC sources in the association catalog with a given SIMBAD 
source type is shown in Table [5j for several different Pid cutoffs. The number of "new" 
associations — RASS/BSC-2MASS/PSC associations not present in SIMBAD — is also given 
for each source type. The new associations include 1464 high quality matches, 679 medium 
quality matches, and 582 low quality matches with previously classified 2MASS/PSC sources. 



J http: //simbad.u-strasbg. f r/simbad/ 
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Of these new associations, only 39 high, 23 medium, and 66 low quality matches are 
with RASS/BSC sources that are classified in SIMBAD. Thus, almost all the new associa- 
tions with previously classified 2MASS/PSC sources provide probable new classifications for 
RASS/BSC sources. For example, the high quality matches provide likely new RASS/BSC 
classifications including 125 (127 if ignoring previous RASS/BSC classifications) double or 
multiple stars, 173 (184) high proper motion stars, and 758 (772) stars. 

Also given in Table[5]is the fraction of 2MASS/PSC sources of a given type that are high 
quality matches. This fraction illustrates the efficiency with which the method described in 
§[2] associates high quality counterparts for a given source type (limited, of course, by any 
selection biases in SIMBAD). Stars, and other bright near-infrared point sources, are among 
the NIR counterparts with the largest fraction of high quality matches, while none of the 
718 quasars or 306 Seyfert type galaxies have > 0.98. This difference is most likely due 
to the definition of LR in the present work, which selects bright NIR sources. 



4.1. Color— Color Diagram 

A RASS/BSC-2MASS/PSC match has flux information both in an X-ray band and 
in multiple NIR bands, and can be plotted with X-ray- J color against J-K color. Such a 
color-color diagram is useful in studying the properties of the matches in the catalog. 

The ROSAT PSPC count rate (0.1-2.4 keV) of each RASS/BSC source is converted to 
an X-ray flux (0.4-2.4 keV) using the HEASOFT tool pimms, by assuming that the spectrum 
is a power-law of photon index 2. Two X-ray fluxes are computed: the uncorrected X-ray 
flux F x (setting A H = in pimms), and the absorption-corrected X-ray flux (Fx)o- The 
value of Ah used to calculate the absorption-corrected flux of each sou rce is derived from th e 



measurements of the Leiden/Argentine/Bonn Survey of Galactic HI (IKalberla et al.l 120051 ). 
using the HEASOFT tool nh. The output of nh chosen in this case is the weighted average 
of the data points from the survey closest to the direction of the RASS/BSC source, and 
it should be noted that this corresponds to the total integrated column to the edge of the 
galaxy. Thus, the applied correction is systematically too large for sources within the Galaxy. 

To calculate the extinction-corrected NIR magnitudes (mj) and (m K )o, it is necessary 
to relate the extinction in the two respective bands to the value of Ah- Such a connection can 
be made by assuming that four separate rela tions between extinction corrections in various 



bands hold for all sources: The first relation (IPredehl fc Schmittlll995l ). 



A H = 1.79xl0 21 A v cm- 2 , (7) 



allows the V-band magnitude correction Ay to be found from the value of Ah; the second 
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relation flSchildlll977f ). 

Ay = 3.2 -E(b-v) > (8) 

allows the B-V c olor correction P(b-v) to be found from Ay; and the final two relations 
jntzpatricklll99fll ). 

Aj = 0.829 P (B -v) (9) 

and 

A K = 0.355 P (B -v), (10) 

allow the J- and K-band magnitude corrections Aj and A^ to be found from Emy). In 
the relation for A-&, the midpoint of the range (0.33-0.38) given for the prefactor is used. 
Having computed the J-band magnitude extinction correction, the corresponding extinction 
corrected flux is 

(Fj) = 5.082xlO- 7 xlO- a4 ( mj - Aj ) , (11) 



wher e the prefactor is the 'In-Band' zero-magnitude flux for the 2MASS J-band (ICohen et al. 



2003|). 



The X-ray to J-band flux ratio against J-K color of the Py > 0.5 matches is shown in 
Fig. [TT] (uncorrected) and Fig. [12] (corrected for X-ray absorption and NIR extinction). The 
source type listed in SIMBAD for the 2MASS/PSC source is shown for source types with 
>50 such P ic j > 0.5 matches, while the remaining classified sources are labeled as 'Other'. 
Sources not listed in SIMBAD, or listed as 'X-ray source' or 'Infra-Red source', are considered 
'Unclassified'. The majority of classified sources are either galaxies or coronally active stars, 
where 'galaxies' include the source types 'Quasar' and 'Seyfert 1 Galaxy', while 'coronally 
active stars' include the remaining source types listed in the joint legend of Fig. [TT]and Fig. [T2] 
(excluding 'Other' and 'Unclassified'). Galaxies and coronally active stars appear to occupy 
distinct regions of the color-color diagram: The galaxies have mj — > (mj — mx)o > 
0.6 and (Fx/ Pj)o > Px/Pj > 3xl0~ 2 , while almost all the coronally active stars have 
(mj- m K )o < rnj — rriK < 11 and Fx/ Fj < (Fx/Fj) < 3xl0~ 2 . There are 3415 unclassified 
Pid > 0.5 matches in the region occupied by coronally active stars, 1647 in the region occupied 
by galaxies, and 442 in the remainder of the color-color diagram. 

The X-ray flux of the bulk of the RASS/BSC sources only spans a range of about one 
order of magnitude, so the large range in Px/Pj seen in Fig. [TT]is mainly due to the large 
range of J-band fluxes of sources in the 2MASS/PSC. Thus, the distribution of matches in 
Fx/ Fj is closely related to the distribution of matches in mj, which can be seen in Fig. [TUJ 
It appears that the brighter population of sources seen in Fig. [10] consists mainly of coronally 
active stars. 
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To determine whether the upper limit on Fx/Fj for coronally active stars is physical or 
due to selection effects (the chosen LR favors bright J-band sources), a search for correlation 
between Fx and Fj was performed for the high quality matches for which the 2MASS/PSC 
source is of one of the stellar source types on the joint legend of Fig. [TT] and Fig. [T5] (i.e. 
not 'Unclassified', 'Other', 'Quasar', or 'Seyfert 1 Galaxy'). Only high quality matches were 
considered to minimize the number of spurious matches, and to exaggerate any selection 
effects. A clear correlation was seen between Fx and Fj along the line corresponding to the 
upper limit on F x /Fj, and the probability th at Fx and Ft are uncorrelated was found to 



be 10 (Spearman Rank-Order Correlation; iPress et al.lll992l ). This correlation indicates 



a physical origin of the upper limit, as no direct dependence on F x is included in LR, but a 
correlation between Fx and Fj could have been introduced through their mutual dependence 
on A^h- However, the probability that Fx and A^h are uncorrelated is 6xl0 -2 , and the 
probability that Ft and Ah are uncorrelated is 0.13 (Spearman Rank-Order Correlation; 



Press et al.lll992l ). so it is unlikely that a correlation between Fx and Fj could have been 
introduced through Ah- The upper limit on Fx/Fj thus appears to be physical in origin, and 
is most likely related to th e saturation at L x / F bol « 10~ 3 seen in rotation-activity diagrams 



for coronally active stars (jGiidell l2004j ) . 



Fig. [13] shows the differential and cumulative number distributions of high quality 
matches that are coronally active stars against Fx/Fj. The sharp cutoff at large Fx/Fj 
is due to the previously discussed physical limit, while the scarcity of very bright J-band 
sources likely affects the shape of the distribution at small Fx/Fj. Since the X-ray flux of 
the bulk of the RASS/BSC only spans a range of about one order of magnitude, the cutoff 
due to the physical upper limit on Fx/Fj for coronally active stars is the most likely cause 
of the similar cutoff seen in Fig. [TOl around mj = 10. Since the 2MASS/PSC is practically 
complete for mj < 15.8 (cf. all the coronally active stars in the RASS/BSC should have 
counterparts in the 2MASS/PSC. The association catalog presented in §3.11 is thus only 
confusion limited for J-band counterparts to coronally active stars, not flux limited. 

The definition adopted for LR favors bright J-band sources (cf. §2.2p . The association 
catalog therefore includes many high quality matches with bright J-band sources, like coro- 
nally active stars, but few matches with fainter NIR sources (cf. Fig. [TUI) . If high quality 
matches with fainter NIR sources are needed, the cross-association can be repeated with a 
different definition of LR. For example, if high quality matches with galaxies are needed, 
LR could be defined to favor NIR objects with (mj — m^jo > 0.6, in much the same way 
bright NIR objects were favored in the present cross-association. 



-20- 



5. COMPARISON WITH OTHER PUBLISHED CROSS-ASSOCIATIONS 

The cross-association method described in §2] is an adaptation to the 2MASS/PSC of 
the method employed in previous work (R00) for the USNO-A2 catalog, so the results of 
the present work can be directly compared to those of the RASS/BSC-USNO-A2 cross- 
association. The quality Q of the USNO-A2 association catalog is 0.652 (R00), so it ap- 
pears that the the fraction of the RASS/BSC for which the present method can identify a 
2MASS/PSC counterpart is is 7±2% larger than the fraction of the RASS/BSC for which the 
previous work could identify a USNO-A2 counterpart. The number of counterparts found in 
both cross-associations is shown in Table [T] for three different Pd cut-offs. The present work 
finds 42% more high quality matches than the RASS/BSC-USNO-A2 cross-association, but 
finds 9% fewer Pd > 0.5 matches. This difference in distribution of Pd values may be due 
to the use of a different passband (NIR) to search for counterparts, which may be sensitive 
to different source classes than the optical passband. The difference may also be due to the 
changes made to LR; in particular, the added p dependence in LR is expected to lead to 
higher pd values for matches in regions with low NIR source density. 

Unlike the RASS/BSC-USNO-A2 cross-association, the present work does not consider 
binary matches. A RASS/BSC source associated with a pair of 2MASS/PSC sources will 
typically have P no -id ~ 0, and appear in the catalog as one of the P id < 0.5 matches. 

Of the high quality matches, 1418 are associations that are not listed in the SIMBAD 
database, and that involve a RASS/BSC source for which a high quality match with a 
USNO-A2 optical source was not presented in previous work. Excluding associations listed 
in SIMBAD and those involving a RASS/BSC source for which a high, medium, or low 
quality USNO-A2 match has been previously been presented, the high quality matches in 
the catalog include entirely new associations for 113 RASS/BSC sources, the medium quality 
matches for 113 sources, and the low quality matches for 808 sources. 

Several cross- associations with other catalogs are described in the RASS/BSC paper 
(jVoges et al. and these are compared to Pd based cross-associations in the RASS/BSC- 



USNO-A2 cross-association (R00). In summary, accounting for source confusion, and not 
selecting counterparts purely based on proximity, make Pd based catalogs more useful than 
the previous cross-associations for classifying RASS/BSC sources. 

Several of the cross-associations listed in Table [T] rel y on plausibility arg uments. One 



such example is a cross association with optical sources (IZickgraf et al.ll2003l ). where each 
optical source in the RASS/BSC error circle was classified through spectroscopy, and the 
"most likely" counterpart then assigned to the X-ray source. This approach does not provide 
a Pd for each match, and relies on previous knowledge of the properties of X-ray emitting 
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source classes. It is therefore unlikely to identify new classes of X-ray sources. 



Machine learning algorithms can be employed to make the plausibility m ethod quantita- 
tive. One application of machine learning algorithms (IMcGlynn et al.ll2004l ) used a training 
sample of previously classified X-ray sources to compute the probability that each RASS/BSC 
source belongs to a given source class. Machine learning algorithms remove the need to spec- 
ify the properties of a given class — the algorithms infer class properties from the training 
sample — and the probability that the RASS/BSC source is of a certain class can be found. 
However, only previously known classes can be used, and the number of objects of a given 
class in the training sample limits the ability of the algorithm to find similar objects. 



6. DISCUSSION AND CONCLUSIONS 

The association catalog presented in §3. II lists each RASS/BSC source, the correspond- 
ing P no _ id , the most likely 2MASS/PSC counterpart, and the P id of the RASS/BSC-2MASS/PSC 
match. The catalog includes 3853 high quality matches (Pid > 0.98), of which < 39 are ex- 
pected to be chance associations with background objects. Of the high quality matches, 
only 1921 are with a RASS/BSC source for which a Pd > 0.98 USNO-A2 counterpart has 
previously been identified (R00). Thus, between the present work and the RASS/BSC- 
USNO-A2 cross-association, there are 4637 RASS/BSC sources with a Pd > 0.98 optical 
and/or NIR counterpart. The present work thus represents a significant improvement over 
previous cross-associations. 

The purpose of uniquely associating RASS/BSC and 2MASS/PSC sources is to obtain 
J-, H-, and K-band fluxes for, and better localization of, RASS/BSC sources, to aid in 
classification. Certain source classes appear to be restricted to specific ranges of X-ray- 
J color and J-K color (cf. §4.ip . so the X-ray, J-, and K-band fluxes of matches can be 
used to select interesting objects for precise classification through follow-up observations. 
Previous classification of the proposed 2MASS/PSC counterpart in a match can also provide 
RASS/BSC classifications. The present work presents likely new classifications for 1425 
previously unclassified RASS/BSC sources, based on the classification of the pd > 0.98 
2MASS/PSC counterpart. 

All the coronally active stars in the RASS/BSC should have counterparts in the 2MASS/PSC 
(cf. §4.ip . so unique association of these RASS/BSC sources with NIR counterparts is con- 
fusion limited. For any coronally active stars in the RASS/BSC that do not presently have 
a known 2MASS/PSC counterpart, better X-ray localizations should allow the 2MASS/PSC 
counterpart to be found. Further, since the fraction of the RASS/BSC sources that do not 
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have counterparts in the 2MASS/PSC is 12 ±1% (cf. §2HJ, there are 2200 ±200 RASS/BSC 
sources with no NIR counterparts in the 2MASS/PSC catalog. These RASS/BSC sources 
are likely associated with NIR sources fainter than the 2MASS/PSC detection limit, and 
represent a sizeable discovery space for classes of X-ray sources other than coronally active 
stars. 

RASS/BSC sources with P no -id ~ 1 have no likely 2MASS/PSC counterparts. Some 
of these sources may be associated with an extended near infra-red source not listed in 
the 2MASS/PSC, some may be associated with NIR sources too far from the RASS/BSC 
position to be identified, and some may be NIR faint objects. The number of RASS/BSC 
sources with P no -id > 0.98 is 334, and these provide a sample which can be searched for 
examples of NIR faint X-ray sources like isolated neutron stars. 
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Fig. 1. — [Left Axis] Differential number distributions of NIR sources as a function of sep- 
aration r between the NIR source and the field center. There is an excess of NIR sources 
in the source fields compared to the background fields. [Right Axis] At each r, the cumula- 
tive distributions from r to 75" are compared using a two-sample K-S test, which gives the 
probability that the same underlying distribution gave rise to both cumulative distributions. 
The K-S probability is less than 10~ 3 for r < 25", demonstrating that the excess in source 
fields is statistically significant (cf. §2.11) . 
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Fig. 2. — [Bottom panel] Differential number distributions in J-band magnitude mj of NIR 
sources in the source and background fields. [Top panel] The difference between the source 
and background field distributions. There is an excess of NIR sources in source fields for 
mj < 17.6, which is statistically most significant for mj < 10.6. 
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Fig. 3. — Normalized number distribution of RASS/BSC sources with A^nir NIR sources in 
the corresponding source field, for iV NIR < 50, and the best-fit model distribution (largely 
obscured by the data distribution). The distribution of RASS/BSC sources has a sharp 
peak around a mode of ./Vnir = 3, and has a tail extending to -/Vnir = 373. There are 1199 
RASS/BSC sources with 50 < iV N i R < 373, with less than 30 RASS/BSC sources at any 
given ./Vnir, which are not included when fitting. The best-fit curve (reduced x 2 = 0.86) 
shown is for 12 ± 1% of RASS/BSC sources being associated with zero NIR sources, 66 ± 3% 
with one, 18 ± 3% with two, and 4.0 ± 0.6% with seven NIR sources (cf. §2.11) . 
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Fig. 4.— [Bottom panel] Number distributions in LR of NIR sources in source 
and background fields. The normalization of the number distributions is such that 
f(N/ field) d(\n(LR)) approximately equals the total number of NIR sources per field for 
that field type (the equality is not exact because of the averaging over bins described in 
$L2}. [Top panel] The calculated R(LR). For \n{LR) < -34.07, R(LR) is set to zero. 
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Fig. 5. — Differential number distributions in R of NIR sources in the source and control 
fields. For R > 0.10 the number density in source fields is above that in control fields, while 
for R < 0.10 the number density in source fields is below that in control fields in all but two 
bins. In total, there is a 9% excess of NIR sources in control fields compared to source fields 
for R < 0.10, and a 78% excess in source fields compared to control fields for R > 0.10 (cf. 
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Fig. 6. — Differential number distributions in Pd for NIR sources in the source and control 
fields. For P i( j > 0.28 the number density in source fields is above that in control fields, while 
for P;d < 0.28 the number density in source fields is below that in control fields in all but two 
bins. In total, there is a 6% excess of NIR sources in control fields compared to source fields 
for P;d < 0.28, and a 222% excess in source fields compared to control fields for pd > 0.28 
(cf. 
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Fig. 7. — [Bottom panel] Normalized cumulative distributions in 1 — P no -id of control fields 
in six different regions of the sky (cf. Table [2]) • See §3] for further discussion. [Top panel] 
Normalized cumulative distributions in 1 — P no -id of RASS/BSC sources in the same six 
regions as for the control fields. The distributions for the regions in the galactic plane follow 
each other closely for 1 — P no -id < 0.75, and lie entirely below the distributions for regions 
outside the galactic plane. 
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Fig. 8. — Normalized number distributions in crp;d of X-ray-NIR matches in the association 
catalog. The high and medium quality matches all have cr Pid < 0.01, while the low quality 
matches have opid values clustered around 0.01, with a tail in the distribution extending to 
0.03. For P i( j < 0.5 almost all matches have opid values clustered around 0.01. See §3. II for 
further discussion. 
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Fig. 9. — Normalized number distributions in separations r of the X-ray-NIR matches in the 
association catalog. The distribution for high quality matches is the most sharply peaked 
(with a mode at 6"), and a broadening is seen with decreasing P^. See §3.11 for discussion. 
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Fig. 10. — Normalized number distributions in J-band magnitude mj of X-ray-NIR matches 
in the association catalog. Almost all high and medium quality matches have mj < 12. In 
the low quality and remaining matches, a population of fainter NIR sources is present. See 
§3. II for discussion. 
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Fig. 11. — X-ray to J-band flux ratio against J-K color of the > 0.5 matches, where the 
point type reflects the source type listed in SIMBAD for the 2MASS/PSC source (cf. legend). 
Only source types for which there are at least 50 members are given a distinct point type, 
while the remaining classified sources are labeled as 'Other'. Sources not listed in SIMBAD, 
or listed as 'X-ray source' or 'Infra-Red source', are labeled as 'Unclassified'. The source 
types are ordered by decreasing number of members (cf. Table ED, and the legend is split 
between this figure and Fig. [12] for legibility (the largest number of sources is 'Unclassified'). 
Two unclassified sources are omitted from the plot because they have mj — itik < —1.5 or 
Fx I 'Fi > 100, and eight other sources are omitted because they do not have an ttt-k listed in 
the 2MASS/PSC. Galaxies and coronally active stars appear to occupy distinct regions of 
the color-color diagram (cf. §4.11) . 
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Fig. 12. — Same as Fig. [TTJ but corrected for X-ray absorption and NIR extinction assuming 
the entire column of material to the edge of the galaxy lies between the earth and the source 
(cf. §4.ip . This correction is expected to approximate the actual correction for extra-galactic 
sources, and is expected to be too large a correction for galactic sources. Six unclassified 
sources with (raj — wk)o < —1.5 or (Fx/Fj) > 100 are omitted, in addition to the sources 
omitted in Fig. [TTJ 
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Fig. 13. — [Right Axis] Normalized number distribution of the 2465 high quality matches 
that are coronally active stars against X-ray to J-band flux ratio. See §4.11 for discussion. 
[Left Axis] Cumulative number distribution of the same matches as for the right axis. 



Table 1. Published catalogs of counterparts for N > 500 RASS/BSC sources 



IXclcI CI1L.C 




l^r'ACC oooA^iotion c o t" o 1 t~\ct'^ 
KjL (JhS-dhSULldllUIl LcLLdlUJi 


AT (XT,. — O/Xb 
-< v l-' v bkg — /a) 


f TlrOnQ T 1 i T ' QCCA^IOTinTI 

vylllljllel 1UI dbbOClelLlUIl 


Huensch et al. 


(1998) 


Bright Star Catalogue 


580 (< 13 = 2.2%) 


< 90" from X-ray source 


Bade et al. 


(19983 


RASS/BSC, |b| > 20° 


3847 (n/a) 


Plausibility 






Hamburg Quasar Survey 






Voees et al. 


(1999) 


HST-GSC 


9759 (1358 = 13.92%) 


< 24" from X-ray source 


Hiinsch et al. 


(1999) 


Third Catalogue of Nearby Stars 


783 (< 16 = 2%) 


< 90" from X-ray source 


Bauer et al. 


(2000) 


RASS/BSC (> 0.1 counts s- 1 ), 


1556 (51 = 3.3%) 


fno-id < 0.5 






|6| > 15°, J2000 S > -40° 










NVSS 






Makarov & Urban 


(2000) 


Tycho-2 


6954 (n/a) 


< 3c7x from X-ray source 


Rutledee et al. 


(2000) 


USNO A-2, B objects 


2705 (18 = 0.7%) 


P\a > 98% 


Rutledee et al. 


(2000) 


USNO A-2, B objects 


5492 (155 = 2.8%) 


Pid > 90% 


Rutledee et al. 


(2000) 


USNO A-2, B objects 


11301 (2034 = 18%) 


Pid > 50% 


Schwoue et al. 


(2000) 


RASS/BSC (> 0.2 counts s" 1 ), |b| > 30° 


2012 (n/a) 


Plausibility 






NED and SIMBAD 










Optical CCD and spectroscopy 






Zimmermann et al. 


(2001) 


Catalogue of Principal Galaxies 


547 (< 66 = 12%) 


< 100" from X-ray source 










Visual Inspection 


Zickgraf ct al. 


(2003) 


RASS/BSC, \b\ > 30° 


4380 (n/a) 


Plausibility 






Hamburg Quasar Survey 






McGlvnn et al. 


(20041 


USNO-B 


6915 (< 1383 = 20%) 


< 30" from X-ray source 






NVSS and SUMSS 




Machine learning classification 


Torres et al. 


(20061 


RASS/BSC Southern Hemisphere 


1791 (n/a) 


< 2.6<7x from X-ray source 






Hipparcos and Tycho-2 




B- V > 0.6 


Present work 


2MASS/PSC, J-band 


3853 (39 = 1.0%) 


P id > 98% 


Present work 


2MASS/PSC, J-band 


6133 (140 = 2.3%) 


P id > 90% 


Present work 


2MASS/PSC, J-band 


10286 (1583 = 15.4%) 


P id > 50% 



a The catalog that was searched for counterparts, and any restriction on which RASS/BSC sources were considered. 

b The number of candiate counterparts presented, and the expected background contamination in number and as a percentage 
of the number of candidate counterparts. 
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Table 2. Number of source and control fields in each of six different regions of the sky (cf. 





$3]), and the 


average p for the fields in 


each rej 


$ion. 






Region 


Galactic Latitude 


Galactic Longitude 


-^CON 


(Pcon) 


^SRC 


(PSRC) 


Entire Sky 


\b\ < 90° 


0° < / < 360° 


41192 


3.17 


18806 


2.85 


1 


50° < \b\ 


0° < / < 360° 


9168 


0.44 


5194 


0.44 


2 


20° < |6| < 50° 


0° < / < 360° 


17270 


1.06 


7767 


1.06 


3 


\b\ < 20° 


|Z — 180°| < 30° 


2455 


3.27 


1049 


3.00 


4 


\b\ < 20° 


30° <\l- 180°| < 90° 


4910 


4.32 


2036 


3.98 


5 


\b\ < 20° 


90° <\l- 180°| < 150° 


4932 


9.12 


1731 


8.90 


6 


\b\ < 20° 


150° <\l- 180° | 


2457 


13.93 


1029 


14.22 



Table 3. Association catalog 



1RXS 


1RXS SIMBAD Type 


PSPC 


C/S (ff) 


Pno~id 


2MASS 


2MASS SIMBAD Type 


mj 


Pid 


J000000.0-392902 


X-ray source 


0.13 


(0.03) 


0.654 


J00000425-3929005 




11.5 


0.134 


J000007.0+081653 


Galaxy 


0.19 


(0.02) 


0.324 


J00000702+0816453 




13.9 


0.670 


J000010.0-633543 


Seyfert 1 Galaxy 


0.19 


(0.03) 


0.449 


J00000981-6335383 




16.7 


0.352 


J000011. 9+052318 


Seyfert 1 Galaxy 


0.26 


(0.03) 


0.266 


J00001172+0523175 




14.4 


0.680 


J000012. 6+014621 


X-ray source 


0.08 


(0.02) 


0.008 


J00001217+0146173 


Double or multiple star 


8.3 


0.991 


J000013. 5+575628 


Star 


0.12 


(0.02) 


0.029 


J00001346+5756387 




9.8 


0.960 


J000019. 5-261032 


X-ray source 


0.12 


(0.02) 


0.329 


J00001916-2610283 




15.5 


0.347 


J000035. 5-280553 


X-ray source 


0.07 


(0.02) 


0.395 


J00003484-2805459 




12.9 


0.474 


J000038.4+794037 


Star 


0.10 


(0.01) 


0.009 


J00004121+7940398 


Star 


8.9 


0.966 


J000042. 5+621034 


Star 


0.16 


(0.02) 


0.006 


J00004167+6210331 


Star 


6.1 


0.989 



Note. — The association catalog is available in machine readable format in the electronic edition of The Astrophysical Journal. A portion 
of the catalog is shown here for guidance regarding its form and content. 
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Table 4. Number of 2MASS/PSC sources in the association catalog that are classified in 
SIMBAD and/or listed as associated with the RASS/BSC source in SIMBAD 



Pid 


Total 


Associated 
in SIMBAD? 


Total 


Classified a 


Unclassified 13 


> 0.98 


3853 


Yes 
No 


1204 
2649 


1186 
1464 


18 
1185 


> 0.9 


6133 


Yes 
No 


1579 
4554 


1538 
2143 


41 
2411 


> 0.5 


10286 


Yes 
No 


2122 
8164 


2057 
2725 


65 
5439 


> 


18568 


Yes 


2609 


2500 


109 


No 


15959 


3386 


12573 



a Sources with a listed type other than 'X-ray source' or 'Infra-Red 
source'. 

b Sources not listed in SIMBAD, and sources listed as 'X-ray source' 
or 'Infra-Red source'. 



Table 5. Number of 2MASS/PSC sources in the association catalog listed in SIMBAD with a given source type 



Source type in SIMBAD 


N(P ld 


> 0.98) 


N(P ld > 0.9) 


N(P ld > 0.5) 


N(P id > 0) 


AT(P id >0.98) b 

N T otal(Pid>0) 




Total 


New a 


Total New a 


Total New a 


Total New a 


Total New a 



Active Galaxy Nucleus 














7 


5 


46 


41 


0.00 


0.00 


Be Star 


4 


2 


5 


3 


5 


3 


5 


3 


0.80 


0.40 


BL Lac - type object 








o 





2 





11 





0.00 


0.00 


Blazar 


o 





o 











1 





0.00 


0.00 


Blue object 


o 





o 


o 


2 


1 


4 


3 


0.00 


0.00 


Brown Dwarf (M < 0.08Mq) 


o 





o 


o 


1 





1 





0.00 


0.00 


Carbon Star 


o 





o 


o 


1 





1 





0.00 


0.00 


Cataclysmic Var. AVI Her type 


o 


o 


o 


o 


12 


o 


22 


o 


0.00 


0.00 


Cataclysmic Var. DQ Her type 


1 


o 


1 


o 


8 


o 


9 


o 


0.11 


0.00 


Cataclysmic Variable Star 


o 


o 


1 


o 


3 


o 


17 


7 


0.00 


0.00 


Ccphcid variable Star 


6 


o 


9 


1 


9 


1 


10 


2 


0.60 


0.00 


Classical Cephcid (delta Cep type) 


1 


1 


1 


1 


1 


1 


1 


1 


1.00 


1.00 


Double or multiple star 


177 


127 


219 


159 


242 


177 


252 


184 


0.70 


0.50 


Dwarf Nova 





o 


3 





30 


4 


45 


12 


0.00 


0.00 


Eclipsing binary 


7 


2 


7 


2 


9 


3 


9 


3 


0.78 


0.22 


Eclipsing Binary Candidate 


1 





1 


o 


1 





1 





1.00 


0.00 


Eclipsing binary of Algol type 


37 


18 


43 


20 


51 


26 


51 


26 


0.73 


0.35 


Eclipsing binary of beta Lyr type 


27 


17 


33 


22 


37 


25 


37 


25 


0.73 


0.46 


Eclipsing binary of W UMa type 


30 


4 


40 


6 


45 


8 


45 


8 


0.67 


0.09 


Ellipsoidal variable Star 


4 


3 


6 


4 


6 


4 


6 


4 


0.67 


0.50 


Emission-line galaxy 




















3 


2 


0.00 


0.00 


Emission-line Star 


15 


6 


21 


9 


23 


10 


23 


10 


0.65 


0.26 


Flare Star 


99 


28 


147 


52 


168 


63 


173 


64 


0.57 


0.16 


Galaxy in Cluster of Galaxies 




















1 


1 


0.00 


0.00 


Globular Cluster 














1 





2 





0.00 


0.00 


Hcrbig-Haro Object 


1 


1 


1 


1 


1 


1 


1 


1 


1.00 


1.00 


High Mass X-ray Binary 


3 





9 





15 





16 





0.19 


0.00 


High proper-motion Star 


286 


184 


392 


264 


477 


337 


493 


353 


0.58 


0.37 


Low Mass X-ray Binary 


1 





1 





2 





5 





0.20 


0.00 


Low-mass star (M < IMq) 








1 





4 


3 


10 


9 


0.00 


0.00 


Nova 








1 





2 





2 





0.00 


0.00 


Nova-like Star 


2 





3 





6 





10 


2 


0.20 


0.00 


Object of unknown nature 




















2 


2 


0.00 


0.00 



Table 5 — Continued 



Source type in SIMBAD 


N(P ld > 0.98) 


N(P id 


> 0.9) 


N(P id 


> 0.5) 


N(P id 


>0) 


iV(P, d >0.98) b 

N T atal(P,d>0) 




Total New a 


Total 


New a 


Total 


New a 


Total 


New a 


Total Ncw a 



Planetary Nebula 








1 





1 





1 








00 


0.00 


Possible Quasar 




















6 








00 


0.00 


Pre-main sequence Star (optically detected) 


15 


1 


44 


3 


64 


10 


71 


10 





.21 


0.01 


Pulsating variable Star 


9 


3 


11 


4 


11 


4 


12 


5 





.75 


0.25 


Quasar 








1 





161 


65 


718 


401 





00 


0.00 


Rotationally variable Star 


12 


4 


12 


4 


13 


4 


14 


4 





.86 


0.29 


Semi-regular pulsating Star 


15 


8 


15 


8 


16 


8 


16 


8 





94 


0.50 


Seyfert 1 Galaxy 








1 





135 


10 


304 


61 





00 


0.00 


Seyfert Galaxy 














1 





2 








00 


0.00 


Spectroscopic binary 


163 


86 


185 


100 


205 


111 


207 


113 





79 


0.42 


Star 


1099 


772 


1586 


1190 


1875 


1439 


2024 


1580 





54 


0.38 


Star in Association 


2 


1 


7 


3 


10 


5 


12 


7 





17 


0.08 


Star in Cluster 


41 


18 


60 


29 


65 


30 


71 


35 





58 


0.25 


Star in double system 


67 


44 


104 


77 


208 


155 


222 


168 





.30 


0.20 


Star in Nebula 














1 


1 


2 


2 





00 


0.00 


SuperNova 














1 


1 


1 


1 


0. 


.00 


0.00 


SuperNova Remnant 




















1 


1 





.00 


0.00 


Symbiotic Star 


4 


3 


5 


4 


7 


4 


7 


4 





57 


0.43 


T Tau-type Star 


94 


2 


185 


12 


238 


15 


247 


17 


0. 


.38 


0.01 


UV-emission source 




















2 


2 





00 


0.00 


Variable of BY Dra type 


79 


35 


91 


41 


107 


49 


109 


51 





.72 


0.32 


Variable of RS CVn type 


110 


16 


123 


19 


136 


24 


136 


24 





.81 


0.12 


Variable Star 


198 


61 


249 


80 


278 


87 


283 


90 





.70 


0.22 


Variable Star of alpha2 CVn type 


2 





2 





2 





3 


1 





67 


0.00 


Variable Star of beta Cep type 


5 


2 


7 


3 


7 


3 


7 


3 





.71 


0.29 


Variable Star of delta Set type 


7 


5 


7 


5 


8 


6 


8 


6 


0. 


88 


0.62 


Variable Star of gamma Dor type 


6 


2 


6 


2 


7 


3 


7 


3 


0. 


.86 


0.29 


Variable Star of irregular type 


5 


5 


6 


5 


6 


5 


6 


5 


0. 


.83 


0.83 


Variable Star of Mira Cet type 


1 





2 


1 


2 


1 


2 


1 


0. 


.50 


0.00 


Variable Star of Orion Type 


3 


2 


8 


4 


14 


7 


19 


9 





.16 


0.11 


Variable Star of RR Lyr type 








2 


1 


3 


1 


4 


2 





00 


0.00 


Variable Star of W Vir type 


1 





1 





1 





1 





1 


00 


0.00 


Variable Star with rapid variations 


1 





1 





1 





1 





1 


.00 


0.00 


White Dwarf 


4 


1 


5 


1 


18 


2 


31 


5 





.13 


0.03 



Table 5 — Continued 



Source type in SIMBAD N(P id > 0.98) N{P ld > 0.9) N(P id > 0.5) N{P ld > 0) 





Total 


New a 


Total 


New a 


Total 


New a 


Total 


New a 


Total 


Ncw a 


Wolf-Rayet Star 


3 





5 


1 


5 


1 


5 


1 


0.60 


0.00 


X-ray Binary 


1 





1 





1 





4 





0.25 


0.00 


Young Stellar Object 


1 





3 


2 


3 


2 


4 


2 


0.25 


0.00 


Young Stellar Object Candidate 




















1 


1 


0.00 


0.00 


All Classes 


2650 


1464 


3681 


2143 


4782 


2725 


5886 


3386 


0.45 


0.25 



a 2MASS/PSC counterparts that are not associated with the RASS/BSC source in SIMBAD. 
b Fraction of 2MASS/PSC counterparts of a given class that are high quality matches. 



